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Abstract 

We calculate the radiative corrections to the current-current correlator of the hy- 
brid current gq{x)"f y iG a ^ v T a q{x). Based on this new result we use the QCD sum rule 
approach to estimate lower bounds on the masses of the J PC =\-+ and 0++ light 
hybrids. 



1 Introduction 

Mesons with exotic quantum number have attracted a great deal of attention in low energy 
strong interaction physics. They constitute another intrinsic construction of matter beyond 
the quark model in QCD. Although these mesons have not been confirmed yet, recent ex- 
periments indeed give some evidence for their possible existence. The E852 Collaboration 
at BNL has reported a J PC = 1 _+ isovector resonance p(1405) in the reaction 7r~p — > r]7r°n, 
with a mass of 1370 ± 16ifg MeV and a width of 385 ± 40t^ 5 MeV 0. This state appears 
to have been confirmed by the Crystal Barrel Collaboration in pp annihilation with a mass 
of 1400 ± 20 ± 20 MeV and a width of 310 ± 50±|g MeV §. E852 lays claim to another 
jpc _ i sovec tor state p(1600) in the reaction -rr~p — > tt + 7t~7t~p, with a mass and width 
of 1593 ± 8 MeV and 168 ± 20 MeV, resp., which decays into pit ||. We will have to wait 
for further confirmation of these states. 



The mass values for the exotic 1~ + state reported by the different experiments disagree 
with most theoretical predictions. The flux-tube model predicts the lowest-lying 1 _+ hybrid 
meson to have a mass of 1.9 GeV 0, which is consistent with lattice QCD studies which 
predict the lightest exotic hybrid 1 _+ to have a mass of 2.0 GeV ||. Besides, the flux-tube 
model also predicts that the h^n and fi7t modes dominate in 1 h hybrid meson decay ||, in 
contradiction with the experiments in which the mode pn is dominant. Differing from the 
flux-tube model predictions calculations based on the QCD sum rule (QCDSR) approach 
seem to be closer to experimental results |7|]. The preliminary results of QCDSR show that 
the lightest exotic hybrid meson has a mass around 1.6 GeV and a dominant pit decay mode 
||. This is consistent with the second 1~ + state claimed by E852. The discrepancies between 
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the different model predictions may come from different sources. The flux-tube model uses a 
non-relativistic linear potential model, which is not so suitable for the light-mass system. On 
the other hand, many effects, such as higher order terms in the OPE and radiative correc- 
tions, may affect the QCDSR predictions. Further theoretical studies are obviously necessary. 

In this paper, we calculate the radiative corrections to the current-current correlator of 
the hybrid current q(x) r y u igG'j lu T a q(x) . Using this new result, we recalculate the masses of 
the 1~ + and ++ hybrids via the standard QCDSR method. We find that, including the 
radiative corrections, QCDSR have less room to accommodate the recent experimental data. 



2 Renormalization of the current = q^f u igG^ u T a q 

The operator-mixing problems associated with the renormalization of composite operators 
have been discussed sometimes ago ||. For instance, a given gauge invariant operator can 
mix with other gauge invariant operators, with non-gauge invariant operators which vanish 
by equations of motion and with operators containing ghosts. The mixing operators must 
have the same CP quantum number and the same dimension as the original one. In our 
case, the complete set of operators which can mix with q^ v igG a T a q is given by 
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Note that there is no dimension-five operator containing ghost fields in the set ([]]). In ([[]) we 
have defined a a p = §[7a, 7/3] and the covariant derivatives = + igA^, = — igA^, 
which act on right and left fields, respectively. The fields and couplings in ([!]) are bare. 
Only and are physical currents. The other currents correspond to so-called nuisance 
operators which vanish by the equations of motion. 

Renormalizing the composite operator corresponding to 3 1 we obtain 

[#] = Zrfl + Z 2 jI + Zsjl + Z^ + Z 5 jI + Z,j% (2) 

In order to determine the coefficients Zi(i = 1,6), we insert the currents into 1PI diagrams 
and extract the ultraviolet divergences. We find that the divergences associated with the 
qqg-vertex are sufficient to determine all counterterms. The relevant diagrams are shown in 
Fig.l. In Feynman gauge, we obtain 

z 1 + 113 92 1 

Zl ~ 1 + T8"16^? 
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z 2 = 

z 3 = 
z, = 

z 5 = 
z 6 = 

where we use dimensional regularization. 
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(3) 
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Next-to leading order correction to the current-current 
correlator 



Let us consider the current-current correlator 

IV(g 2 ) = if*xe*(0\Tfa(x)J+(p)}\0) 



(4) 
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where j M (x) = u(x)j u igG® v T a u(x) . The invariants IL^Q 2 ) and U S (Q 2 ) correspond to the 
contributions from 1~ + and ++ states and their excited states, respectively. The leading 
order contribution to @) including the quark and gluon condensate contributions has already 



been given in ]TT 
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+ 
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(5) 



M=f-> 



There is no difference for isovector and isoscalar currents in this order. The next-to- 
leading order correction to the perturbative part of n„(g 2 ) and n s (g 2 ) can be obtained by 
calculating the Feynman diagrams in Fig. 2, where Figs. 2m and 2n only contribute to the 
isoscalar states. The technique of the calculation which we use here was firstly proposed in 
0- 

Let us now briefly comment on the radiative corrections to the correlator of the current 
9Qlhlv'iG < t L ,T a q with and l + ~ quantum numbers. The isovector current correlator has 
the same radiative correction as the isovector current correlator of the (0 ++ , 1 h ) current. 
However, the isoscalar current correlator has different radiative corrections since diagrams 
Figs. 2m and 2n, which correspond to the mixing with pure gluonic states, now give zero 
contributions. 
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Each of the diagrams Figs.2i-2n is gauge-parameter independent by itself because the 
current is antisymmetric in the Lorentz indices. We have checked on gauge invariance for 
the sum of diagrams Figs.2a-2h by doing the calculation in a general covariant gauge and 
found that the result is gauge-parameter independent. Explicit Feynman gauge results for 
the diagrams are listed in the appendix. 

In the present application only the isovector currents are of interest. In the MS-scheme, 
the next-to leading corrections to the correlator of isovector currents is given by 
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where we take the light quark mass to be zero for convenience. Although (||) include con- 
densates which are proportional to the light quark mass, their contributions are very small 
compared to those of operators with the same dimension. Throughout this calculation taking 
zero quark mass is a good approximation. 

Eq.(|6|) does not contain the complete next-to leading order correction, one also has to in- 
clude the contribution from the renormalization of the current. By inserting the renormalized 
currents (0) into the correlator 



d 4 xe^2Z«°(Q\T{jl(x),ji + (0)}\0), (7) 



we obtain 



IW) = -^(g 2 ) 3 ln(^)((-- + -n / )^ + (---n / )^ln(^)) (8) 
vKqj 240vr 3W ; V 2 16 40 " n V 36 6 }J tx V 2 

4807T' ;i fi z 144 120 7i 36 6 tt fi z 

In Eq.(0) we have introduced the notation Z^ s = ZiZ g — 1 and Z" a = Z^ii = 2, ..6). For 
the original hybrid current y one also has to include the counterterm Z g resulting from the 
renormalization of the bare coupling constant g. The final result for the radiative corrections 
to the correlator is finally obtained by adding the contributions of (0) and (Bp. 



4 Sum rules for 1 + and ++ hybrid mesons 

For sum rule applications we need the spectral density associated with the current-current 
correlator. The spectral density p v (s) = ImU v (s) is defined via the standard dispersion 
relation 

iyg 2 ) = ^ r°ds J v(s) 2 , + j:a k ( q r, (9) 

tx Jo s n (s - q 2 ) ^ 
where the a k are appropriate subtraction constants to render Eq.@ finite. 
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From (^j), ([|) and (||) we obtain 

. . ot s (u) o,„ 1301 q s (u) 17a s (u) , , s 1 .. „ , 

= ^ 1+ w^-^V )+ i s((asG)+8fts(mMU))(1 ^ 

m 1 480vr 2 v 2160 7T 36 tt V v ; 

-(y (mim) + - )a - ^-(a s G 2 ) - ^m(MG«), 

where we have set n/ = 3. 

On the other hand, the spectral density is saturated by narrow physical resonances and 
the continuum. We therefore write 



q ^Y- - 9fiu)Pv{s) = ^(0|j M |i?)(i?|>|0)7r5(s - m 2 R ) + continuum, (12) 

1 R 



where we have assumed that the mass tur is much larger than the decay width of the hybrid, 
so that the imaginary part of the propagator has been replaced by tt5(s — m 2 R ). 

In order to extract information on the lowest-lying resonance, it is usually assumed that 
the lowest-lying resonance dominates the spectral density. The contribution of higher excited 
states can be suppressed by applying the Borel transformation Lm to both sides of Eq. @. 
One thus has 

Ro(t) = ML M Tl v (q 2 ) = - / e~ s ! M p v (s)ds. (13) 

7T JO 

The upper limit of the intergral can be replaced by a finite number sq. The contributions 
beyond the threshold so are considered to result from the continuum. Rq(t) is the zeroth 
moment. Higher order moments are defined by R k = (M 2 ^) k R (M). Resonance masses 

9 Ft 

can be obtained by taking the ratio m R = with the assumption that only a single 
narrow resonance dominates. In principle any value of k can be chosen to determine the 
resonance masses. However, since we have to truncate the series of the power expansion, 
using higher moments will damage the convergence of the OPE. Besides, it is also arbitrary 
to define the scalar function Il k u (q 2 ) = j^Ii^q 2 ) by extracting a general tensor factor 

{q 2 ) k {^-9^ (14) 

from the correlator Although ( |TB"D can be considered as a higher moment of U k (q 2 ), to 
the order of OPE that we are considering, the sum rules for Il k (q 2 ) and U v (q 2 ) are obviously 
different. For instance, as pointed out in [11], the dimension-six operators of (||) do not 



contribute to sum rule (0), while they play an important role in stabilizing the sum rule of 
njj(g 2 ) in the case of k = 1. We shall consider the two cases k = and k = 1 in turn. 
The single particle matrix elements contributing to ([3~2]) are parametrized as 



(0|^|V) = ie,f v ml (15) 

In the narrow resonance approximation the sum rules are independent of the matrix element 
(|15|). When the decay width of the resonance is comparable with its mass, the approximation 
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fll^D is no longer valid. Then the information contained in (|TJ|) may become important. We 
will comment on this later on. By using (ITOj)- (p~5|) we obtain 



for the sum rule (13) and 



m, 



Jo°e ST sp V)S (s)ds 
Io° e- ST p v , s (s)ds 

f*° e- ST Pv (s)ds 



(16) 



/ S0 e-^ Pv {s)ds/s - ^a s (uu) 2 - j^g 3 (G 2 ) + ^a s m(uGu) 
f*°e- ST p s (s)ds 



So e- ST p s {s)ds/s - 8ir 2 a s (uuy 



for the sum rule for U^(q 2 ) with k = 1 , where the spectral densities p v ,s{s) are given in (|T0|). 
The various parameters entering in ([T6|) - (pTT|) are specified as ]12| 



K QCD = 0.25GW, 



m = O.OlGeV, 



m(uu) 



t2 2 



(a s G 2 ) = 0.04Gev\ g 3 (G 3 ) = l.lGeV 2 (a s G 2 ), f n = 0.132GW 

a.(^) = — ; „2 r , A* = 2GeV, g(uGu) = 1.5Gev 2 (uu) 

9 In ) 

In Fig. 3 and Fig.4 we show a mass plot of the 1 h state in its dependence on the Borel 
parameter M in the two sum rule ratios. The second sum rule gives a smaller mass which 
we will take as the lower bound. 

The sensitivity of the mass to the choice of the threshold value Sq is obvious. The mass 
will go to infinity when both s and M go to infinity, because, when M goes to infinity, 
the Borel measure ([13]) no longer suppresses the continuum. In order to give a reasonable 
estimate, we set sq around 4 GeV 2 . Fig. 5 shows that m v takes values around 3 GeV when 
So is set to infinity. In Fig. 6 one cannot find any stable point in the two-dimensional (sq,M) 
space. Therefore, before fixing sq, we cannot make any precise prediction for the hybrid 
mass. However, if one believes that the 1 h hybrid mass lies around 2 GeV, sq should 
be larger than 4 GeV 2 . This results in a lower bound for the 1 h mass of 1.55 GeV(see 
Fig. 2). The radiative corrections enhance the lower bound which gives QCDSR less room to 
accommodate recent experimental data. 

The prediction of the mass is also sensitive to the form of the spectral density. This is 
mostly due to a truncated OPE. The contributions from higher dimension operators are very 
likely not small. The uncertainty from the narrow resonance approximation in ([121) does not 
seem to reduce this discrepancy. This can be checked by replacing the narrow resonance 



form tt5(s — m|j) in (12) by the Breit-Wigner form 



m 



2 \2 



+ Y 2 m 2 



(19) 



R 



Tr is the width of the 1~ + hybrid. When we choose Yr = 200MeV and a parametrization of 
the matrix element ([TS]) in the form (0^1 V) = ie^mnf v s x with f v constant and x = 0.5-=-1.5, 
the hybrid mass is not sensitive to using the full Breit-Wigner form propagator. We show 
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the change in the Figs. 3 and 4 by dotted lines. The discrepancy between the two cases k=0 
and k=l is still big and the mass predictions become somewhat larger. 

The sum rule for the ++ hybrid is shown in Figs. 7 and 8 for the two cases, k=0 and k=l 
resp., where we set s = 7 GeV 2 . The radiative corrections reduce the discrepancy between 
the two cases. Similar to the 1 _+ case, when s goes to infinity, the prediction of mass is 
around 3 GeV. It means that the contribution of the continuum dominates over that of the 
resonances. Therefore, the value chosen for the threshold sq is important. 



5 Summary 

In summary, we have calculated the next-to-leading order corrections to the two point corre- 
lator of the current gq r ) v iG a VLV T a q[x). We recalculated the masses of the 1 h and ++ hybrids. 
We find that the radiative corrections reduce the lower bound of the 1 h mass and leave less 
room for QCD sum rules to fit the recent experimental data. 

Note added in proof: While preparing this paper for publication, we became aware of 
a recent paper by K. Chetyrkin and S. Narison which addresses similar problems |13| . 
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Appendix In this appendix we list the results of calculating the diagrams Fig. 2 in the 
Feynman gauge for the correlator (HI). 



Fig.2a 
Fig.2b 
Fig.2c 
Fig.2d 
Fig.2e 
Fig.2f 
Fig.2(g + h + i) 
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Fig.2l 
Fig.2m 
Fig.2n 
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We use the abbreviation C = ^^(q 2 ) 3 ln(^). 
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Figure captions 

Fig. 1 Feynman diagrams for the renormalization of the hybrid current. Dots stand for 
the current vertices. 

Fig. 2. Feynman diagrams for the next-to-leading calculation. Dots stand for the current 
vertices. 

Fig. 3. 1 h hybrid mass m v versus Borel variable M for the first sum rule Eq.fllBl). 
Dashed line gives the result for the leading order calculation. Solid line includes radiative 
corrections. Dotted line gives the results when using the Breit-Wigner resonance propagator. 

Fig. 4. 1~ + hybrid mass m v versus Borel variable M for the second sum rule Eq . (|I"7|) 

Fig. 5. 1~ + hybrid mass m v versus Borel variable M for the second sum rule Eq.([T7|) 
when s goes to infinity. 

Fig. 6. Three-dimensional figure of 1~ + hybrid mass m v vs. the Borel variable M and 
the threshold parameter sq. 

Fig. 7. ++ hybrid mass m s versus Borel variable M for the first sum rule Eq. (O) . Solid 
line inculdes radiative corrections. 

Fig. 8. ++ hybrid mass m s versus Borel variable M for the second sum rule Eq. (|17D . 
Solid line includes radiative correction. 
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Fig.l 
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